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>■ ■ The Taylor-Reynolds and Reynolds number {Rex and Re) dependence of the 

^jC) • dimensionless energy dissipation rate q = eL/uf^^^^ is derived for statistically sta- 

Q I tionary isotropic turbulence, employing the results of a variable range mean field 

theory. Here, e is the energy dissipation rate, L is the (fixed) outer length scale, 
—I. ■ and Ui^rms a rms velocity component. Our fit-parameter free results for c^{Re\) 



CK ■ and also for Re\{Re) are in good agreement with experimental data. Using the 

i?e-dependence of c^ we account for the time dependence of the mean vorticity uj{t) 

for decaying isotropic turbulence, which was recently experimentally examined [M. 

Smith, R. J. Donelly, N. Goldenfeld, and W. F. Vinen, Phys. Rev. Lett. 71, 2583 

j^ ■ (1993)]. The lifetime of decaying turbulence, depending on the initial Rex^, is pre- 

^ ■ dieted and found to saturate at 0.18L^/t/ oc Re\ q {v is the viscosity) for large Rexfl- 

H : PACS: 47.27.-i, 47.27.Gs 

•i-H - 



Dimensional analysis has proven to be a powerful tool in turbulence research, 
giving a number of key features of turbulent spectra |I|, |^. The main idea is to 
connect small scale quantities such as the energy dissipation rate e with the large 
scale quantities such as the outer length scale L and a rms velocity component, 
Ui^rms- More precisely, following Richardson's cascade picture of turbulence p, 0, 
it is argued that for fully developed turbulence 

e = Ceul,,„,jL (1) 

holds, where q is a dimensionless constant in the range of 1. 

However, it is not clear, for what Reynolds number Re the large Re limit (|I|) is 
reached. For small Re the dimensionless number c^ clearly depends on Re. E.g., for 
plane Couette flow with shear 2ui^rms/L in one direction we have Ce{Re) = 4/ Re. 
Here and henceforth, following Sreenivasan [Q , we defined the Reynolds number Re 
as 

Re = Lui^rms/j^, (2) 

where z/ is the viscosity. Even for large Re it is not clear whether the i?e-dependence 
of Ce vanishes. Recent experiments on the Taylor Couette flow, which is bounded, 
seem to show that a logarithmic i?e-dependence persists up to very high Re 0. 
(Here, Re is defined with the velocity of the outer rotating cylinder.) It can be 
accounted for by employing Prandtl's boundary layer theory [§, P, 0, yet math- 
ematically only Ce < const can be proven up to now [§, 0. For unbounded flow 
experiments favor q = const for high Re, see Sreenivasan's collection of data for 
grid turbulence 0]. 

Another way of expressing the physical contents of (0) is to give the Re- 
dependence of the Taylor-Reynolds number Rex = Xui^rms/^, where A = 
ui,rms/{diUi)rms IS the Taylor microscale. With e = 15z/((9iUi)^^^, valid for isotropic 
turbulence 0, and eqs. (|I|), (^) we get 



Rex = Jl5Re/c,{Re). (3) 



If we have q = const for large Re, we thus also have Rex oc v i?e. 

In this Letter we will first derive explicit expressions for c^{Re) and c^{Rex) for 
unbounded flow. They only depend on the Kolmogorov constant b, which is known 
from many experiments to be 6 = 8.4 0]. We do not introduce any free parameter. 
This is possible as we employ the results from the variable range mean field theory 
|[T0|| , which embodies the Navier-Stokes dynamics. We thus offer a way to go far 
beyond dimensional analysis. In the second part of the paper we apply our results 
to decaying turbulence, which has very recently been experimentally examined and 
analyzed by dimensional analysis by Smith, Donelly, Goldenfeld, and Vinen ]TT| . 
Both our results for c^{Rex) and for the time dependence of the mean vorticity u!{t) 
in decaying turbulence are in good agreement with experimental data. 



We start from the final result of Efiinger and Grossmann's variable range mean 
field theory |T^. For homogeneous, isotropic turbulence a differential equation for 
the velocity structure function D{r) = {{u{x + r) — u{x)y) is derived [0], namely 




(4) 

For small r the solution is D{r) = er^/(3z/), whereas for large r we get the well known 
Kolmogorov result D(r) = b^er^^^. The Kolmogorov constant b can be calculated 
within the approach of ref. [|10|. But rather than taking the mean field result b = 6.3 



Tod , we will use eq. (^) with the experimental value b = 8.4 0, as this value gives 



better agreement with the measured structure function D{r), see Fig. 4 of ref. |TD| . 
Eq. (^) is an energy balance equation. The first term in the brackets corresponds to 
viscous dissipation, the second one can be interpreted as eddy viscosity. 

Integrating the differential equation (^) from to the outer length scale L we 
get 

D^L) + 3bhuD{L) - bh^L^ = 0. (5) 

Here we have assumed that (^ holds for all r up to r = L, which is definitely not 
the case for bounded flow and even for unbounded flow slight corrections might be 
necessary, see below. Anyhow, if L is large enough and the flow is isotropic, it is 
D{L) = 2{u^) = 6u1j.^g. With eqs. (|ID and (||) the quadratic (in e) equation (H) 
can be written in dimensionless form as 



18 /6 
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which is easily solved to give 



Ce{Re)=c,,oo<{^4-] -^ + \l^+^-^l^.}- (7) 





The dependence of c^ on the Tay/or-Reynolds number Rex can be obtained from 
eqs. (ID and (H), 
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In both formulae we have introduced q^oo = (6/6)'^'^ = 0.60 only for convenience; 
the Re- and Rex- dependences are purely determined by the Kolmogorov constant 
b. 

For large Re or Rex, the function Ce{Rex) indeed becomes constant, Ce{Rex) = 
Ce,oo = (6/&)^'^^, as Sreenivasan finds for grid turbulence with biplane square mesh 
grids 0. The experimental value q^oo ~ 1.0 is slightly larger than our result c^^oo = 
0.60. The reason for the slight discrepancy is likely due to non universal boundary 
effects. We have assumed (H) to hold for all r up to r = L, thus D{L) = Guf^^^ = 



6(eL)^/^, which aheady results in Cj^oo = (6/6)^/^ = 0.60. Yet for r around L, D{r) 
will not scale as D(r) = 6(er)^/^ and we have D{L) < b{eLy^^. This non universal 
boundary effect might be treated by introducing an effective, geometry depending 
5(e//) < b instead of b in eq. (g), defined by D{L) = 6('=^^)(eL)2/3. (Note that an 
introduction of b^^^^^ already in (|^) is inappropriate, as boundary effects should not 
be seen in D{r) for small r.) To get the experimental 0^,00 = 1-0, one should have 

^(e//) _ g/^2/3 = Q < fo = 8.4. 

The quantity Ce{Rex)/ce^oo, eq. (||), is plotted in Fig. 1, together with Sreeni- 
vasan's experimental data for grid turbulence. For Re\ ~ 50 the function Ce{Re\) 
saturates at c^^oo, in good agreement with the data. 

For really small Re eq. (|^) can - strictly speaking - no longer be applied, as 
laminar flow is never isotropic, whereas eq. (^) only holds for turbulent, isotropic 
flow. Note that Fig. 1, starting with Rex = 5, does not include the laminar case (as 
seen from the inset), since in laminar flow Rex looses its meaning. If we perform the 



small Re hmit Re < ^J3b^/8 = 14.9, Rex < ^Jhh^^/A = 27.2, nevertheless, we can 
get Ce{Re) = 18/ Re, independent oi b, as expected, since b characterizes the highly 
turbulent state. The oc -Re~^-dependence is correct. The prefactor 18 is - again as 
expected - too large, if compared to the highly anisotropic laminar Couette flow 
with shear in only one direction, see above. In more isotropic laminar flow, e.g., in 
flow with shear in three directions, the agreement for small Re will be better. 

Equations (^) and (0) give the function Rex{Re), see inset of Fig. 1, which 
strongly resembles latest experimental measurements for grid turbulence by Castaing 
and Gagne |l^, [l^. Also Grossmann and Lohse obtain a very similar curve Rex{Re) 



from a reduced wave vector set approximation of the Navier-Stokes equations |T 



Here, for large Re we have Rex = Jl5{b/Q)^/'^\/Re ~ byRe. For small Re it is 



R^\ = \/5/QRe. Extrapolating these two limiting cases, the crossover between them 
takes place at Reco = 18(6/6)'^/^ = 29.8, corresponding (via eq. (|^)) to Rex,co = 
21.6, which seems very realistic to us, cf. also ||13|| . 

Up to now we applied our theory to statistically stationary turbulence. But it also 
offers an opportunity to analyze c^ecaymg' turbulence. Most experiments on decaying 
turbulence have been performed in wind tunnels up to now, where the distance from 



the grid gives the decay time t, if the mean velocity is known |T^. In this kind of 
experiment the outer length scale L grows with time t, as the wake behind the grid 
becomes wider with increasing distance. Yet in a very recent new type of experiment 
performed by Smith, Donelly, Goldenfeld, and Donelly |Tl| , L can be kept fixed. In 



that experiment a towed grid generates homogeneous turbulence in a channel filled 
with helium II. The decay of the mean vorticity uj{t) is measured by second sound 
attenuation [^, W\. As in [^ we assume that Navier-Stokes dynamics can be 



applied to this fluid, see the discussion in [11 



In the theoretical analysis of this experiment, for very high ReQ = Re{t = 0), 
Rexfi = Rex(t = 0), the quantity Cg can be considered constant. But the smaller 
Re{t) or Rex{t) become with increasing time t, the more important are the cor- 



rections seen in eqs. (0), (H). The total energy (per volume) of the flow is 
E = uf.^g/2 = 3ulj.^g/2. The decay of the fully developed turbulence is governed 
by the differential equation 



u 



2\3/2E3/2(t) 



'^-' c,(E(t))(^) ^_A-=-E(t). (9) 



The outer length scale L is fixed, as in the experiment we are referring to [TT|]. We 



change variables to i?e(t) = J2/?>LJ E{t)/u and obtain 

Re = --^CeiRe)Re^. (10) 

Integrating eq. (|TD|) with the initial condition Re{t = 0) = Reo gives the time 
dependence Re(t) of the Reynolds number, 

t 3 rR<t) -dx , , 

11 



Here, for simplicity we have introduced the viscous time scale r = 1? jv and the 



constant 7 = wSfo^/S = 9/ce,oo = 14.9. The integral can be calculated analytically. 
We define the indefinite integral as 



^^^'^ = 2^? {-^ + a/TT^R^} + ^ log I 



7 + V7^ + Re' 



(12) 



Re J 

Thus the time dependence of Reif) is given by the inverse function of 

t{Re)lr = 3{F{Re) - F{Re^))/c,^^. (13) 

Note again that there is no free parameter in eq. ([13|) , all quantities on the rhs can 
be expressed in terms of the Kolmogorov constant h. 

Imagine now the limiting case of large Re^ and also large time, but Re{t) < 
Reo still large, i.e., t not too large. For large Re both terms in (|1^) contribute 
(2i?e)-^[l + 0(7/i?e)]. Thus F{Re) = 1/Re and from (0) we exphcitly get 

«,w^(i +£ipi)-«X(i)-. (M) 

Viieo 3 T/ Ce,oo Vr/ 

In Fig. 2 we plotted Re\{t), calculated from eqs. ([T3| ) and (^, for several Re^. The 
scaling law Re\{t) ^ 3v^(t/r)" /q^oo (corresponding to (0)) only starts to be 
observable for Reo ~ 10'^, i.e., Re\fl ~ 156. 

In the final period of decay, i.e., for very large t (large enough so that Re{t) <^ 7), 
we get Re{t) = 27V^exp (-6t/r) and E{t) = l9/A)eb^iy'L-' exp {-Ut/r). An 
exponential decay for very large t also holds for decaying turbulence with growing 
outer length scale L(t) []T5|. 



To compare our results with the hehum II experiment |Tl[], we have to calculate 
the mean vorticity a;(t). Vorticity always causes strain in the flow. It can be shown 
0] that uu^ = e. Thus 

uu^ = e = -E = -3u^ReRe/L'^. (15) 

With eq. {^) we get 

Tto{t) = Re^Rec,{Re) = ^c:^Re{t)sf^+^^+R^), (16) 

where the universal law on the rhs again only depends on the Kolmogorov constant 
h and on the time t. 

Next, we estimate the lifetime ti of the decaying turbulence. For this purpose 
we calculate, how the Kolmogorov length [El 



^{t) = {v^/e{t)f" = .fM^), (17) 

depends on time t. Of course, ri{t) will increase with time, as the turbulence becomes 



weaker and weaker. The behavior can be obtained from eqs. (0) and (13) for any 
Rcq. If Reo is large enough, scaling r]{t)/L oc (t/rY^^ can develop. 

How to define the lifetime ti of the turbulence? As the crossover in the structure 
function D{r) between the viscous subrange and the inertial subrange happens at 
r ^ lOrj [^, we define the lifetime ti by the condition 10ri{ti) = L. A Reynolds 
number Rei is associated with this time ti via the eqs. ([T6| ) and ([T7|) . We calculate 
Rei = 20.3, Re\^i = 16.0, which is, as it should be, near to the viscous-turbulent 
crossover in the curve Rex{Re), which occurs at Reco = 29.8, Re\co = 21.6, 
see inset of Fig. 1. With this definition we obtain the lifetime ti of the decaying 
turbulence for any given Reo (or, via (|), Rex^) as 

tiiReo)/r = 3iF{Rei) - F{Reo))/c,,oo- (18) 

We plotted ti{Rexfl) in the inset of Fig. 2. For small Rexfl the lifetime ti{Rexfl) 
grows logarithmically with Rexfl- For very large Rex,o it saturates at ti{Rex,o) = 
2>TF[Rei)/ce^oo = 0.18r, i.e., it becomes independent on Reo, if measured in time 
units of r = L^/z/. The lifetime, if measured in seconds, of course increases oc Reo oc 
Re\ in the limiting case. 



Finally we compare with the data of the helium II experiment ||rT|. First we have 
to embody the boundary effects, as they should be larger in the helium II experiment 
than in grid turbulence, because the turbulence decays in a tube. Indeed, in ref. 
HH '^i^oo^ = 36.4 is given (using our definition of q, eq. (|I])), corresponding to 

5(e//) 4 Q/(ci%f^)m = 0.55. 

Using this 6(^^-^) instead of b (or ^(^■/'/), c'fl^\ respectively) in eqs. (|12D and (p!6D, 
we plotted uj{t) in Fig. 3, together with Smith et al.'s experimental data |lT|. The 
two curves show the same features. For small t there is no power law. For medium t 



the theory gives Tuj(t) ^ 3^/^(t/r)~^/^/c^'^^-'^). The power law exponent —3/2, which 
is clearly seen in the experimental data, has already been derived by dimensional 
analysis ||ll|. But in theory the — 3/2-power law for uj{t) extends much further 



than observable in experiment, because there experimental noise hinders observation 
already for t > 10s. For further comparison a reduction of the experimental noise is 
essential. Experiments with dramatically increased sensitivity of the detectors are 
in progress [^ . 

The theoretical lifetime ti of the decaying turbulence (calculated with b^^^^^ = 
0.55) is ti = O.OlSr = 140s, much larger than the 10s in which the a;(t)-signal can be 
measured. Thus viscous effects, arising from the i?e-dependence of q for smaller Re, 
see eq. (|^), only become important for a time t/ ^ 10s. So the slight decrease in the 
measured uj(t) signal for t ^ 10s is not due to them, as one might have thought, but 
possible due to the uncoupling of the normal and superfluid components of helium 
II (which was used as the fluid in the experiment |ll|]), as speculated in [p!T| . 



We summarize our main results. We first calculated the functions c^{Re\) and 
Rex{Re), eqs. (^, (|^, and (||), from a variable range mean field theory [^, which 
goes far beyond dimensional analysis. We then applied our results to decaying tur- 
bulence, highlighted by the expression for the time dependence oi uj{t), eq. (|1^). All 
results are in good agreement with experiment. To even improve the agreement, fu- 
ture work has to be done to embody non universal boundary effects in this approach. 
A way to do so is to introduce fe^^'^^) instead of b in eq. (^). Alternatively, one could 
get rid of the boundary effects by calculating a high passed filtered velocity field 
from the experimental data, so that the non universal effects are filtered out. 
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Figures 

Figure 1: The i?eA-dependence of Cf^/Cf^^oo according to eq. (||) with h = 8.4. Also 
shown is Sreenivasan's collection of experimental data |^. The inset shows Re\{Re) 
according to eqs. (|^) and (|^). The crossover from the laminar scaling Re\ oc Re to 
the turbulent range scaling Re\ oc \fRe takes place at Reco = 30, Re\^co = 22. 



Figure 2: Time dependence of the Taylor- Reynolds number Re\{t) for 

freely decaying turbulence with fixed outer length scale for different Rcq = 
W\ 10^ 10^ 10^ {Rexfi = 1-58 ■ 10^ 4.98 ■ 10^ 1.58 ■ 10^ 498), left to 
right, cf. eqs. (|^), (0), and ([131) with b = 8.4. The dashed line denotes the lifetime 
ti i.e., the time, when viscous effects start to dominate, see text. The inset shows 



the i?eA,o-dependence of the lifetime ti of the turbulence, calculated from eqs. (|T8D, 
(|ll), and d) with b = 8.4. 



Figure 3: Time dependence of the mean vorticity uj{t) from Smith, Donelly, Gold- 



enfeld, and Vinen's experiment |Tl|] for -Rca.o = 634, solid line. In that experiment 
L = 1cm, V = 8.97 ■ 10~^cm^/s, so that r = L'^/u = 1.11 ■ lO^s. The dashed line is 
our prediction for uj{t) according to eq. (^) with b^^^^^ = 0.55. The corresponding 
lifetime ti = 0.013r = 140s is given, too. 
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